Force extension curves for polymers with hydrophobic patches: Role of Disorder 
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Motivated by single molecule experiments on biopolymers we explore force vs. extension behavior 
of flexible chains with hydrophobic segments using numerical simulations and theory. We find that 
in addition to the fraction of hydrophobic patches their spatial distribution along the backbone of 
the chain play a major role in altering its mechanical response. These results are discussed in light 
of the helix-coil model for biopolymers. 
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INTRODUCTION 

Biopolymers such as DNA, RNA, and proteins carry 
out a majority of cellular processes that are crucial to sus- 
taining life[l[. These processes are often associated with 
conformational changes of biomolecules. Several single 
molecule experiments, that probe mechanical properties 
of biomolecules, aim to relate primary structure to their 
native conformation and hence function^. 

Theoretical investigations on mechanical properties of 
proteins, DNA and other bio-macromolecules often use 
coarse-grained descriptions modeling structural compo- 
nents as non-linear elastic springs. Sequence hetero- 
geneities are manifested as stiff and compliant elements 
along the chain backbone and are modeled by a site de- 
pendent bending modulus Ki (where i corresponds to the 
lattice site). Statistical mechanical techniques rooted in 
the helix-coil model Q and the worm-like chain model 
H,[f| have been employed to explore stretching and bend- 
ing response of biopolymers These models however 
assume an annealed form of disorder, allowing for the in- 
terconversion between the "helix" and "coil" phases of 
the polymer chain function of either temperature 
or force. How does the mechanical response of a poly- 
mer chain with quenched disorder differ from that of an 
annealed one? 

Motivated by this question we explore the mechanical 
properties of a polymer chain with hydrophobic and hy- 
drophilic beads representing the quenched disorder using 
Langevin dynamics simulations and theory. We inves- 
tigate the stretching response of a polymer subject to 
equal and opposite tensile forces applied to its ends. In 
our simulations the polymer chain is composed of hy- 
drophobic and hydrophilic segments with the ability to 
sclf-hybridise. Our main result is that the mechanical 
response of the polymer chain is determined both by the 
fraction of hydrophobic contacts, /, as well as their ar- 
rangement along the chain backbone. We consider three 
disorder configurations: (a) random, corresponding to 
hydrophobic segments uniformly distributed throughout 
the chain backbone, (b) periodic, where a repeat unit 
of hydrophobic beads are intersperced with hydrophilic 
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FIG. 1: Schematic representation showing hydrophobic 
groups (black circles) arranged along the chain backbone in 
(a) random, (b) periodic, and (c) block copolymer arrange- 
ments respectively. The fraction of hydrophobic beads is 
/ = 0.5. 



ones at regular intervals, and (c) block copolymer, where 
the hydrophobic groups arc clumped together at one end 
of the chain. A schematic representation is shown in 
Figure Q] and the stretching response of different config- 
urations are shown in Figure [3J 

Khokhlov and co-workers have investigated de novo 
protein design of A — B copolymers that give rise to 
folded configurations [1-15|. The goal of this study was 
to identify arrangements of A and B molecules on the 
chain backbone that gives rise to globular states akin to 
native protein conformations. In order to investigate this 
a bead-spring polymer is equilibrated in a poor solvent 
forming a globule. The atoms that form the core shell 
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of this globular structure are "colored" [HJ 1 1 11 and the 
chain is stretched out to note the position of these atoms 
along the chain backbone. A set of backbone atom ar- 
rangements that lead to optimal packing geometries can 
then be identified. However the mechanical response of 
chains with imposed quenched disorder to an applied ten- 
sile force has not been explored. 

Our results can be summarized as follows: we observe 
a sharp first order transition in the chain size at a crit- 
ical force F c for all three disorder realisations at a high 
value of the hydrophobic bead fraction / (see figure [3]). 
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For a fixed hydrophobic fraction / (and / > 0.5) the 
critical force F c for the block copolymer distribution is 
observed to be higher than that of the periodic and ran- 
dom distribution. Interestingly, the difference between 
the critical forces corresponding to periodic and random 
distributions is negligible. The second result is that the 
critical force is zero for all three distributions for low val- 
ues of /. A non-zero critical force F c is observed above 
a threshold value of / which is found to be the same for 
all three distributions. The nature of the transition is 
first order like for block copolymer distribution whereas 
appears to be a crossover phenomena for the random and 
periodic ones. 



SIMULATION METHOD 

The co-polymer is modeled as a freely-jointed bead- 
spring chain consisting of N beads linearly connected by 
N — 1 springs. The co-polymer comprises of two types 
of monomers (or beads), they are either hydrophobic or 
hydrophillic. In stating the extent of hydrophobicity of 
the monomers, we are emphasizing the relevance of our 
coarse-grained simulations to dilute solutions of polymers 
in an aqueous solvent. As illustrated in Fig. [TJ we con- 
sider three different arrangements of hydrophobic beads 
along the polymer chain. The hydrophobic beads are dis- 
tributed either in a random, periodic or block co-polymcr 
arrangements. Figure Q] represents the three different dis- 
tributions for a chain with a fraction, / = 0.5, of hy- 
drophobic groups. 

Each bead in the polymer chain represents a Brownian 
particle, and the length of each link (spring) corresponds 
to the Kuhn length, Iq, for the polymer. The dynamics 
of each bead is described by the Langevin equation, 



dt 2 



(1) 



where mi, and Vj represent the mass, position and ve- 
locity of particle i, respectively. The random force, F£(i), 
arising from the bombardment of the monomer bead by 
the solvent molecules is described by the fluctuation- 
dissipation theorem, 



(Fn*)FK*')> = 6{k B T6(t-t')5 
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where fcs is the Boltzmann constant, T the absolute tem- 
perature, and £ is the frictional drag. Externally applied, 
equal and opposite forces, Ff xt , act on the terminal beads 
of the chain. However, Ff xt = for the remaining poly- 
mer beads. Equation [T] was solved numerically using the 
velocity- Verlet algorithm, with a time step At = O.OOlr, 
where r = lo(m/ e^,/) 1 / 2 . 

The net interaction potential, £/(ry), is given by, 



U(rij) = U LJ + U bond , 



(3) 



corresponding to the excluded volume interactions and 
bond stretching, respectively. The excluded volume in- 
teraction between any two beads is described by, 

tj ( n _ / 4e LJ [(cr/r lJ ) 12 - (a/rij) 6 ], m < 2.5cr, , , 
Ulj V*> ~ \ 0, ry > 2.5<r, [ ) 

where = [r^ — Yj , clj is the Lennard- Jones interac- 
tion parameter and a is the bead diameter. In our sim- 
ulations, we have set the bead diameter a = 0.75/o flij ]. 
We use €lj, lo and m as scales for energy, length, and 
mass, respectively. The mass of a polymer bead is set to 
a value of one. The bond stretching potential between 
adjacent beads is given by the FENE potential, 



1 - (r b /R max y 



(5) 



where r b is the separation distance between adjacent 
beads, R ma x = l-5Zo is the maximum allowable sep- 
aration distance between bonded beads and k is the 
spring constant. The simulations were executed in 
the open-source molecular dynamics simulation package 

LAMMPS[l7|. 

We used a chain size of N = 100 in our simulations. 
The polymer chain is enclosed inside a cubic box of edge 
length 50 lo with periodic boundary conditions imposed 
in all three directions. In order to model the appropri- 
ate effective solvent interactions, we set e^j — 1.5 for 
the hydrophobic beads and e^j = 0.1 for the hydrophilic 
ones. For a homopolymer, a choice of e^j = 0.1 for all 
beads simulates a polymer chain that demonstrates self- 
avoiding walk statistics, whereas a value of e^j = 1.5 
results in a collapsed globule. When e^j = 0.4, we re- 
cover Gaussian statistics for the homopolymer. We var- 
ied the hydrophobic bead fraction, /, between and 1, 
to explore the effects of disorder on the force-extension 
response of the polymer chain. 

The initial configuration of the chain was such that a 
fraction, /, of the beads were assigned as hydrophobic 
and the rest, 1 — /, were considered to be hydrophilic. 
Initially, the polymer chain is equilibrated in solution for 
a period of 4 x 10 6 time steps in the absence of any exter- 
nal force (F ext = 0). The chain was then stretched from 
the equilibrium state to full extension by an incremental 
increase in pulling force, F ext , from a value of 0.1 to 10. 
The pulling cycle is followed by a retraction cycle where 
pext j g ^ eCT easec j f rom 10 to 0.1. During the pulling and 
the retraction cycles, the polymer is equilibrated for a 
period of 5 x 10 6 time steps at a given value of F ext . We 
calculate the time-averaged value of the radius of gyra- 
tion, R g , over the equilibration period. Consequently, the 
extension of the polymer at that value of applied force 
is calculated. The data was averaged over Nrun = 5, 
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FIG. 2: Figure shows the average radius of gyration (R g ) 
over different disorder realisations as a function of the hy- 
drophobic groups for (a) block copolymer (filled circles o), 
(b) periodic distribution (filled triangles A), and random dis- 
tribution (open diamonds 0)- The average is over 5 different 
realisations for the random system. 



disorder realizations for the three types of distributions 
and each value of /. In the case of the random distribu- 
tion, the five realizations correspond to five statistically 
different arrangements of hydrophobic groups along the 
polymer backbone, for the same value of /. The results 
obtained from these different realizations were used to 
calculate the average values of and standards for and cor- 
responding to a particular value of force and fraction of 
hydrophobic groups 



RESULTS AND DISCUSSION 

Figure [2] describes the variation of the radius of gyra- 
tion of the polymer, R g , as a function of /, for all three 
distributions at zero applied force. We find that as / 
increases, R g decreases monotonically for all three cases. 
An increase in / results in aggregation of hydrophobic 
groups leading to a monotonic decrease in R g . The vari- 
ation of R g with / follows near identical trends for both 
the random and periodic distributions for the entire range 
of /. Both distributions display a sigmoidal behavior, 
where the polymer chains show a transition from a coiled 
state to a collapsed, globular state for / > 0.4. This is 
in contrast to the behavior of R g with / for the block 
copolymer distribution. For / < 0.2, R g values for the 
block copolymer distribution are close to the ones corre- 
sponding to the periodic and random ones. However, for 
0-2 / 0.8, the rate of change of R g with respect to / 
is lower than the other two distributions. Also, the block 
copolymer does not go through a sharp phase transition 
as the other two distributions. For / > 0.8, R g values 
for all three distributions tend to a common value. At 
/ = 1.0, R g is the same for all three distributions. 



In the case of a block copolymer, the hydrophobic 
groups at one end form a globule, while the rest of the 
chain assumes a coil structure. As / increases, a greater 
fraction of the chain is in the globular state, whereas the 
hydrophilic block remains as a coil. Consequently, the 
entire chain does not form a compact globule for any 
value of / apart from / = 1. However, for the periodic 
and random distributions, hydrophobic groups are dis- 
tributed along the entire polymer chain, which enables 
attractive interactions between hydrophobic groups that 
are far separated along the chain contour. This results 
in a first order coil-to-globule transition for / « 0.35. 

1 /2 

Figure [3] shows extension, (R 2 ) /L, vs. force, F, 
for all three distributions at four different hydrophobic 
fractions / = 0.2, 0.4, 0.6 and 0.8. For small hydropho- 
bic fraction, / = 0.2 (Fig. [3^a) ) , the chain extension in- 
creases monotonically with F for all three distributions, 
with little variation between them. For higher hydropho- 
bic fractions (/ = 0.4, 0.6, 0.8), a first order globule to 
coil transition is observed at a critical force, F cr . In the 
case of periodic and random distributions, the hydropho- 
bic groups, which are distributed along the entire chain 
contour, hybridize with each other to form a compact 
globule. Therefore, when the magnitude of the applied 
force is small, the globule is only slightly distorted from 
its equilibrium conformation. Consequently, the exten- 
sion remains nearly constant for small F. In contrast, 
the block copolymer undergoes a relatively larger exten- 
sion for smaller values of F (F < F cr ) for / = 0.4, 0.6, 0.8. 
This is because in poor solvents, the block copolymer con- 
sists of a hydrophobic globule at the end of a flexible, hy- 
drophilic coil that is easy to stretch under the application 
of a tensile force. Once the long wavelength fluctuations 
along the chain backbone have been smoothened out, the 
end-to-end distance of the chain remains constant upon 
an increase in the force F. However, beyond the critical 
force F cr , the chain unravels undergoing a jump discon- 
tinuity in size. On further increase in F above F cr , the 
chain stretches out to full extension. As is expected, the 
large force behavior is identical for different disorder re- 
alisations as shown in Fig. [3] 

Figure |4] shows the variation of F cr as a function of 
/ for all three distributions of hydrophobic groups. For 
/ < 0.3, the critical force, F cr , required to uncoil the 
polymer chain is zero for all three distributions. This 
is because for small /, the total number of hydropho- 
bic contacts are few. Thus, the polymer remains in a 
coil-like state, and starts stretching continuously even for 
small values of F. In case of the block copolymer, F cr in- 
creases discontinuously to a large value of approximately 
3.2 (from a zero value) for / rj 0.35, and with increasing 
/, F cr quickly saturates to a value of 4.0. In contrast, 
for the periodic and random distributions, F cr appears 
to increase linearly with / for / > 0.35. It is interesting 
to note that the F cr vs. / curves arc nearly identical 
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FIG. 3: Figure showing force F vs. extension {R g )/L for dif- 
ferent values of the hydrophobic fraction / and different dis- 
order realisations along the chain backbone, block copolymer 
(filled circles o), periodic (triangles A), and random distribu- 
tion (diamonds <))■ 




FIG. 4: Figure shows the variation of critical force F cr as a 
function of the fraction of hydrophobic groups / for differ- 
ent disorder realisations (a) block copolymer (circles o), (b) 
periodic distribution (triangles A), and random distribution 
(diamonds 0). The average is over 5 different realisations for 
the random system. 



for the random and periodic distributions. For / = 1 
the critical force of uncoiling F cr = 4 (in dimensionless 
units) is same for all three distributions. This is because 
all three distributions have the same structure when all 
the beads are hydrophobic. 



SUMMARY 

We have utilized Langevin dynamics simulations to 
investigate the force vs. extension behavior of a sin- 
gle co-polymer composed of hydrophobic and hydrophilic 



groups arranged in periodic, random and block distribu- 
tions along its backbone. Our simulations demonstrate 
that the mechanical response of a heteropolymer under- 
going a force-induced globule-coil transition is dependent 
on the distribution of hydrophobic groups along its back- 
bone. All three disorder realizations form globular struc- 
tures at large hydrophobic fractions, and show a first or- 
der transition in their sizes at a critical force, F c . How- 
ever, for a given fraction, the critical force is observed 
to be higher for the block copolymer in comparison to 
the periodic and random distributions. In addition, our 
simulations show that the force-extension responses of 
the periodic and random distributions are nearly iden- 
tical. From the force-extension curves, we plotted the 
critical force, F c , as a function of the hydrophobic frac- 
tion, /. Since the polymers behave like coils for small 
hydrophobic fractions, F c = for all three distributions. 
An important result is that F c is non-zero for / 0.35 
for all three distributions. On further increase in the hy- 
drophobic fraction, the critical force shows a first order 
transition for the block copolymer. However, it increases 
continuously for the random and periodic distributions. 
Interestingly, we again note that the transition behav- 
iors of the random and periodic distributions are nearly 
identical. 

We now place our results in context of existing the- 
ories on the subject. The collapse transition of disor- 
dered polymers has been the subject of several stud- 
ies [19l - l21| . Statistical techniques rooted in replica meth- 
ods has been used to explore globular conformations 22 1 
(and references therein). Further optimal distributions 
of hydrophobic and hydrophilic groups along the primary 
sequence leading to folded conformations have been in- 
vestigated in context of rational design of proteins 0- 
15l . Il8| . and statistical mechanical theories of random 
systems 21 1 have been employed to explore such con- 
formational states of A-B/H-P heteropolymers. 

However mechanical behavior of disordered polymers 
have been less investigated [23J] . Our results illustrate the 
difference in mechanical properties of such heteropoly- 
mers having similar bulk structural properties but dif- 
fering in the local arrangement of disorder undergoing 
a force induced globule-coil transition. Transfer matrix 
calculations based on the helix-coil model for annealed 
disorder 0, 0] can be adapted to the case of quenched 
disorder. However this is omitted from the present paper 
for brevity and the fact that it does not yield to any new 
insights over those obtained from simulations. 

We hope that this work will spark interest amongst 
polymer chemists interested in rational design of poly- 
mers that aim to identify sequences that lead to specific 
folded structures. Further, we believe this work would 
also be useful to single molecule biophysicists interested 
in understanding protein conformational transitions and 
allostery. 
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